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Abstract: Climate warming has and will continue to exert a significantly impact on crop phenology in the 
past and coming decades. Combining observed data of phenology and a crop growth model provides a 
good approach for quantitating the effects of climate warming and crop management on crop phenology. 
The purpose of this study is to determine the impacts of temperature change, sowing date (SD) 
adjustment and cultivar shift on maize phenology in northern China during 1981—2010. Results indicated 
that climate warming caused anthesis date (AD) and maturity date (MD) of maize to advance by 0.2—5.5 
and 0.6-11.1 d/10a, respectively. Due to climate-driven changes in maize phenology, three growth periods 
of maize, i.e., vegetative growth period (VGP; from sowing to anthesis), reproductive growth period (RGP; 
from anthesis to maturity) and whole growth period (WGP; from sowing to maturity) shortened by 0.2-5.5, 
0.4-5.6 and 0.6-11.1 d/10a, respectively. With SD adjustment (ic, SD advancement), AD and MD 
occurred early by 0.5-2.6 and 0.1-3.4 d/10a, respectively. SD adjustment caused duration of VGP of 
maize to prolong. However, duration of RGP slightly shortened by 0.1-1.3 d/10a. Furthermore, due to 
cultivar shift, MD of maize significantly delayed by 4.9-12.2 d/10a. Durations of VGP, RGP and WGP of 
maize prolonged by 0.2—4.1, 1.6-8.4 and 4.3-11.8 d/10a, respectively. In conclusion, our results indicated 
that cultivar shift, to some extent, could mitigate the negative impact of climate warming on maize 
phenology. 
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1 Introduction 


Long-term observation and related climate study have demonstrated that the earth's climate has 
undergone significant changes (IPCC, 2013). Climate warming has and will continue to exert a 
major impact on agricultural production (Piao et al., 2010; Pongratz et al., 2012), resulting in a 
food security problem on the regional or global scale (Lobell et al., 2008; Tao et al., 2012). 
Phenology is the study of the timing of recurring biological events, as affected by biotic and 
abiotic factors, such as climate warming (Ma et al., 2012). Crop phenology, an indicator of crop 
responses to climate variability and change, is mainly driven by the prevailing climate conditions 
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during its growth season (Xiao et al., 2013; Wang et al., 2015; Liu et al., 2017). Therefore, 
variations in phenological periods of a crop are strong indicators of changes in environment and 
climate (Estrella et al., 2007). Generally, climate warming accelerates the growth and 
reproduction of crops and eventually has a negative impact on crops' production (Xiao and Tao, 
2014; Tao et al., 2014a). In addition, crop phenology is also affected by field management 
practices, such as cultivar shift and sowing date (SD) adjustment (Xiao et al., 2016a, b; Rezaei et 
al., 2017; Liu et al., 2018). Therefore, cultivar shift and SD adjustment could reduce the negative 
effect of climate warming on crop production (Xiao et al., 2013; He et al., 2015; Liu et al., 2018). 

In China, maize produces raw material for an array of multiple products and is one of the most 
important crops (Tao et al., 2014b). The total output of maize grain increased 12-fold during the 
past six decades. The average grain yield of maize ranks the second among cereal crops (Li, 
2009). A number of studies documented that an increase in climate warming during growth season 
had a negative impact on maize phenology during the past several decades (Li et al., 2014; Tao et 
al., 2014a; Xiao et al., 2016a; Liu et al., 2017; Abbas et al., 2017; Guo et al., 2018). However, the 
negative impact of climate warming on maize phenology can be mitigated by crop management 
(Liu et al., 2013; Wang et al., 2016; Xiao et al., 2016a). For example, developing and introducing 
new cultivar with a longer duration of growth period had a positive impact on maize phenology 
and ultimately enhanced yield under climate warming (Sacks and Kucharik, 2011; Xiao et al., 
2016a). Therefore, to deeply understand the processes and mechanisms controlling maize 
responses to ongoing climate warming, researchers must investigate the changes of maize 
phenology (Liu et al., 2013; Xiao et al., 2016a). 

Recently, impacts of climate warming on crop phenology are of great concern and have been 
extensively evaluated by various methods, such as experimental observations, statistical models 
and crop simulated model (Tao et al., 2014a, b; Xiao et al., 2016a, b; Liu et al., 2018). In general, 
due to the complex relationships between climate, cultivar shift and management, the effects of 
climate change and improvement of crop management on crop phenology are difficult to be 
clarified by field experiments and/or statistical models (Wang et al., 2016). However, by using 
crop simulated model, researchers can easily separate the impact of multiple and interacting 
factors on crop phenology into single factor impact (Liu et al., 2013; Xiao et al., 2016a; Wang et 
al., 2016). Therefore, crop simulated models are widely used to quantify the effects of climate 
warming and crop management on crop phenology (Holzworth et al., 2014). For example, 
agricultural production systems simulator (APSIM) is a software system which allows models of 
crop and pasture production, residue decomposition, soil water and nutrient flow and erosion to 
be readily re-configured to simulate various production systems and soil and crop managements 
to be dynamically simulated using conditional rules (McCown et al., 1996). 

In this study, observed data from 19 agricultural meteorological stations in northern China (NC) 
during 1981-2010 were used to study the response of maize phenology to climate warming. 
Furthermore, the APSIM model was used to separate the effects of climate warming, SD 
adjustment and cultivar shift on maize phenology. Thus, the main purposes of this study are (1) to 
investigate the variation in maize phenology in NC during 1981-2010; (2) to correlate the 
observed and stimulated phenological period with climate warming; and (3) to quantify the effects 
of temperature change, cultivar shift and SD adjustment on maize phenology. 


2 Materials and methods 


2.1 Materials 


Data of maize phenology were collected from 19 agrometeorological stations (here referred as to 
19 stations) in NC during 1981-2010. The stations, operated by the China Meteorological 
Administration (CMA) and provincial meteorological administrations, are distributed across the 
main maize cultivated regions in NC (Fig. 1). These stations distributed in Hebei (HB), Shanxi 
(SX), Gansu (GS) provinces and Xinjiang Uygur Autonomous Region (XJ) and Inner Mongolia 
Autonomous Region (IMG) (Fig. 1). Phenological parameters including dates of sowing, 


emergence, heading, anthesis and maturity of maize were recorded during growth season. 

Furthermore, a total of 5 agrometeorological stations (Fengning station in HB, Changzhi station 
in SX, Wuwei station in GS, Tacheng station in XJ and Chifeng station in IMG; here referred as to 
5 stations) that had detailed records about cultivars during 1981—2010 were selected to investigate 
the effects of climate warming and SD management on maize phenology in NC (Fig. 1). Old 
maize cultivars were replaced by improved cultivars by local farming community within 3-5 
years (about 9 cultivars per station during 1981-2010) (Table 1). The long-term detailed 
phonological records during 1981—2010 provide us an opportunity to look insight into the 
response of maize phenology to SD adjustment, climate warming and cultivar shift (Xiao et al., 
2016a). 
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Fig. 1 Locations of 19 investigated agrometeorological stations in northern China. IMG, Inner Mongolia 
Autonomous Region; HB, Hebei Province; SX, Shanxi Province; GS, Gansu Province; XJ, Xinjiang Uygur 
Autonomous Region. Five stations with underlines were selected for analyzing the cultivar shift. 


Table 1 Maize cultivars grown at 5 stations in northern China during 1981-2010 


Number Station Cultivar 
1 Chifen Jidan-101, Huangmo-417, Zhongdan-2, Beiyu-3, Yedan-4, Nongda-3138, Shuoqiu-8, 
E Nongda-364, Danyu-202 and Zhengdan-958 

2 Fenenin Jidan-101, Feng-71, Zhongdan-2, Jinghuang-417, Fengdan-3, Tangkang-5, 
anme Nongda-3138, Nongda-95 and Denghai-20 

3 Changzhi Nanban-12, Zhongdan-2, Nongda-60, Yudan-13, Nongda-3188, Tunyu-7, Sanbei-6, 

Dafeng-5 and Jiyu-811 
4 Wuwei Zhangye-488, Jindan-2, Zhongdan-2, Liangdan-1, Shendan and Shuoqiu-9 
5 Tacheng Hongmaya-3, Shuangjiao-42, Maya-2, Jinhuanghou, Xinyu-3, Yidan-7, Yidan-8, 


SC-704 and Denghai-3672 


Moreover, climate data including mean (T), minimum (Tmin) and maximum (Tmax) temperatures, 
sunshine duration and precipitation at 5 stations during 1981-2010 were used. For each station, 
daily solar radiation (Ra) was calculated from sunshine duration at 5 stations using the 
Angstrom—Precott (A—P) equation (Prescott, 1940): 
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n 
Ruzar z) R, >’ (1) 


where Ra is the extra-terrestrial solar radiation (MJ/(m?.d)); a and b are the AP coefficients; and n 
and N are the actual and theoretical sunshine durations (h), respectively. In this study, a and b 
coefficients for 5 stations were calibrated using observed radiation value. 


2.2 Methods 


We calculated mean temperature (T) of three growth periods, i.e., vegetative growth period (VPG; 
from sowing to anthesis), reproductive growth period (RGP; from anthesis to maturity) and whole 
growth period (WGP; from sowing to maturity) based on observed dates of sowing, anthesis and 
maturity at each station. Moreover, we calculated accumulated temperature (AT, the sum of daily 
temperature >0°C) of VGP, RGP and WGP based on observed phenological periods. By using 
regression analysis, we calculated the trends of maize phenology, durations of growth periods and 
climate variables (T and AT) during 1981—2010. The trend of each variable was determined by 
using the following linear regression model: 


Y, =kX, +b > (2) 
where Y; is the observed phenological date or climate variables (T or AT) in the year i; k is the 
linear regression slope; b is the slope intercept; and X; is the year i (i=1, 2, 3, ..., 30). Bivariate 


correlation analysis was performed to determine the correlation between duration of growth 
periods (VGP, RGP and WGP) and T and AT. Then, statistical significance was determined with 
two-tailed t-test analysis. 

Generally, observed phenological periods (Mov) were mainly impacted by climate warming, 
cultivar shift and SD adjustment. Using the APSIM model, we set two experiments (Mı and M2) 
to determine the effects of climate warming, cultivar shift and SD adjustments on maize 
phenology at 5 stations. In experiment Mı, SD was set as the average value of observed SD during 
1981-2010. The design of Mı experiment is favor to detect the trend in phenological periods 
affected by climate warming (Table 2). In experiment M2, SD was set as the observed SD during 
1981-2010. The design of M2 experiment is favor to determine the trend in phenolgical periods 
affected by climate change and SD adjustment (Table 2). The impacts of cultivar shift and SD 
adjustment on maize phenology were determined via linear regression analysis as in Equation 2. 
Here, Y=Movi—M2; (for the impact of cultivar shift) and Y=M2;-Mi; (for the impact of SD 
adjustment), where Movi, Mi; and Mo; are the values of each variable in the year i, respectively. 


Table 2 Details of experiments used to determine the effects of climate warming, cultivar shift and sowing date 
(SD) adjustment on maize phenology in northern China 


Experimental 


cade Explanation Effect 
Mo Observed SD during 1981-2010 Climate warming, cultivar shift and SD adjustment 
Mı Average value of observed SD during 1981-2010 Climate warming 
M: Observed SD during 1981-2010 Climate warming and SD adjustment 
M:—M, Difference between M> and M, SD adjustment 
Moø-M2 Difference between observed SD and M, Cultivar shift 


2.3 APSIM model 


In this study, APSIM version 7.7 was used to separate the effects of climate change and SD 
management on maize phenology. A detailed description of APSIM model and its crop modules 
are recorded in Holzworth et al. (2014). APSIM model requires daily climate data as an input, 
including solar radiation, maximum and minimum air temperatures and precipitation. APSIM 
model was calibrated and validated by using observed data from 5 stations. Firstly, the cultivar at 
each station during 1981—1985 was identified. Then, observed data during 1981—1983 were used 
to calibrate the parameters of APSIM model. Also, observed data during 1984—1985 were used to 
validate the model. Finally, validated model was run on climate data during 1981—2010 to 
simulate station-based dates of anthesis and maturity. 

Overall, the simulated dates by APSIM model and observed dates in the field of anthesis and 


maturity fitted well at 5 stations (Fig. 2). The differences between simulated and observed dates 
of anthesis and maturity were less than 5 d, suggesting that APSIM model performed well in 
simulating maize phenology in this study. 


Fig. 2 Fitting curve of simulated dates by APSIM (agricultural production systems simulator) model and 
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observed dates in the field of anthesis (a) and maturity (b) at 5 stations in northern China. DOY, day of year. 


3 Results 
3.1 Variation in maize phenological period during 1981-2010 


In NC, maize was generally sown from early April in GS province to early May in IMG region 
(Table 3). Anthesis generally occurred from mid-July (DOY (day of year) 191) to early August 
(DOY 214), and maturity was generally occurred from late August (DOY 239) to late September 
(DOY 272) in NC (Table 3). As shown in Table 3, the durations of VGP, RGP and WGP were 


Table 3 Observed sowing date (SD), anthesis date (AD) and maturity date (MD) and durations of three growth 


periods (VGP, RGP and WGP) of maize at 19 stations in northern China during 1981-2010 


Province/region Station SD (DOY) AD (DOY) MD (DOY) VGP (d) RGP (d) WGP (d) 
IMG Chifeng 119+6.1 206+4.0 259+6.9 8745.0 5347.1 140+7.9 
Naiman Banner 1186.0 2063.8 260+10.0 8746.6 548.6 142+12.3 

Tongliao 124+4.7 209+4.0 26242.9 8544.6 5343.5 1385.1 

Tuquan 127+6.2 21445.5 26745.2 8747.6 544.7 1407.9 

Zalantun 134+4.6 2113.8 2645.1 7144.7 536.0 13145.8 

Jungar Banner 120+6.0 205+4.7 265+4.3 8545.9 605.2 145+7.0 

HB Fengning 1135.1 2066.0 2627.7 936.5 5644.4 149+8.5 
Huailai 1174.1 202+4.4 2592.8 85+4.2 563.4 1413.8 

SX Xiyang 124+5.8 2115.9 272+4.3 86+6.0 6144.2 148+4.8 
Xixian 120+7.2 2104.2 2676.8 907.0 5747.0 1478.1 

Xinzhou 114+4.0 202+2.7 2574+4.9 8844.5 55+4.4 14345.3 

Changzhi 118+4.4 2053.6 262+4.4 8745.2 5743.5 144+5.9 

GS Chengxian 100+45.2 191+4.2 23947.2 9145.3 4845.5 139+46.3 
Jiuquan 10743.8 20343.9 266+6.2 974.5 637.2 159+7.5 

Tianshui 108+4.9 2016.9 255+7.7 9345.7 543.6 146+6.6 

Wuwei 101+4.0 2105.2 2686.1 109+6.2 58+4.8 167+6.4 

XJ Tacheng 12248.7 200+7.5 2567.2 7846.9 565.6 1347.3 
Usu 1135.8 196+5.0 25245.4 8244.6 563.9 1395.0 

Xinyuan 115+7.6 2007.3 2546.2 85+6.5 544.7 13947.2 

Average 117+48.5 2055.7 2607.4 88+7.0 5643.4 1448.3 


Note: IMG, Inner Mongolia Autonomous Region; HB, Hebei Province; SX, Shanxi Province; GS, Gansu Province; XJ, Xinjiang Uygur 
Autonomous Region; VGP, vegetative growth period; RGP, reproductive growth period; WGP, whole growth period; DOY, day of year. 
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78-109, 48-63 and 131-167 d, respectively (Table 3). This result indicated spatial differences in 
the lengths of growth periods at 19 stations due to the diverse climate conditions. 

As shown in Figure 3a, SD of maize during 1981-2010 advanced at 17 of 19 stations, with a 
significant advancement at 6 stations (P<0.05). AD of maize occurred early at 14 of 19 stations, 
and the trend was significant at 5 stations (P<0.05; Fig. 3b). As shown in Figure 3c, MD 
advanced at 10 of 19 stations, and the trend was significant at 1 station (P<0.05). MD delayed in 
9 stations (mainly in ING region and HB province), and the delay was significant at 4 stations 
(P<0.05; Fig. 3c). 
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Fig. 3 Variations in sowing date (a), anthesis date (b), maturity date (c) and durations of three growth periods 
(vegetative growth period, reproductive growth period and whole growth period) (d—f) of maize in northern China 
during 1981-2010. Black circles show that the trend is significant at 5% probability level. 


The advance or delay in maize phenology induced the changes in durations of growth periods. 


As shown in Figure 3d, VGP shortened in 9 stations (mainly in XJ and GS) and the trend was 
significant at 2 stations (P<0.05). In contrast, VGP lengthened in ING and SX, and the trend was 
significant at 5 stations (P<0.05; Fig. 3d). Also, RGP lengthened at 16 of 19 stations and the trend 
was significant at 6 stations (P<0.05; Fig. 3e). As shown in Figure 3f, WGP lengthened at most of 
the stations and the trend was significant at 7 stations (P<0.05). 

At 19 stations, SD advanced by 1.8 d/10a, AD slightly advanced by 1.2 d/10a, and MD slightly 
delayed by 0.5 d/10a (Fig. 4). As shown in Figure 4, the durations of VGP, RGP and WGP 
lengthened by 0.6, 1.8 and 2.4 d/10a, respectively. 
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Fig. 4 Trends of maize phenology under different phenological periods at 19 stations in northern China during 
1981-2010. SD, sowing date; AD, anthesis date; MD, maturity date; VGP, vegetative growth period; RGP, 
reproductive growth period; WGP, whole growth period. Bar indicates standard error. The lower and upper box 
boundaries indicate the 25" and 75" percentiles across 19 stations, respectively; whiskers below and above the 
box indicate the 10" and 90" percentiles across 19 stations, respectively; and the line and small box within each 
box indicate median and mean values, respectively. 


3.2 Variations in T and AT during 1981-2010 


At 5 stations, variation in T during three growth periods of maize showed an increasing trend 
(Table 4). Variation of Tvcp was significant at Fengning and Wuwei stations, variation of Trop 
was significant at Wuwei and Tacheng stations and variation of Twop was significant at Wuwei 
station (Table 4). Furthermore, variation in AT during ATvcp, ATrcp and ATwo of maize 
significantly increased (P<0.05) at 5 stations. These results implied the climate warming trends 
during growth periods of maize. And the shift of maize cultivars preferred to the species with a 
higher AT requirement for coping with climate warming. 


Table 4 Variations in temperature (T) and accumulated temperature (AT) during three growth periods (VGP, 
RGP and WGP) of maize at 5 stations in northern China during 1981—2010 


Station Tver (°C/a) Trop (°C/a) Twap (°C/a) AT var (°C/a) ATrep (°C/a) ATwap (°C/a) 
Chifeng 0.011 0.012 0.011 4.32T° 8.907" 13.277 
Fengning 0.038" 0.006 0.023 7.496" 4.717" 12.206™ 
Changzhi 0.001 0.006 0.003 6.156" 0.248 6.374" 
Wuwei 0.066” 0.074" 0.070% 4.748" 8.571" 13.303” 
Tacheng 0.005 0.077" 0.039 1.291 15.402”° 16.697" 
Note: * and * denote significant differences at P<0.05 and P<0.01 levels, respectively; VGP, vegetative growth period; RGP, 


reproductive growth period; WGP, whole growth period. Tygp means the temperature during vegetative growth period. The above 
abbreviations in Table 4 are the same combination of T or AT and different growth periods. The abbreviations are the same in Table 5. 
As shown in Table 5, the duration of three growth periods was negatively correlated with T. 
The results suggested that a warming temperature could reduce phenological period of maize. 
However, the duration of three growth periods was significantly and positively correlated with AT 
(Table 5), indicating that cultivar shift with a higher AT requirements had significant impacts on 
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maize phenology. 


Table 5 Correlation between duration of three growth periods (VGP, RGP and WGP) and temperature (T) and 
accumulated temperature (AT) at 5 stations in northern China during 1981-2010 


Station Tver Trap Twop AT var AT Rap AT wep 


Chifeng —0.274 -0.518"* —0.295 0.801" 0.941” 0.852" 
Fengning 0.181 0.312 —0.351 0.830" 0.722" 0.808" 
Changzhi —0.054 -0.478" —0.508™ 0.875” 0.677” 0.764" 

Wuwei —0.368 —0.324 -0.421° 0.744" 0.694" 0.519" 
Tacheng -0.515™ —0.165 0.251 0.894" 0.819%" 0.735” 


Note: * and ** denote significant differences at P<0.05 and P<0.01 levels, respectively. VGP, vegetative growth period; RGP, 
reproductive growth period; WGP, whole growth period. 


3.3 Effects of climate warming, crop management and cultivar shift on maize phenology 


At 5 stations, SD occurred early by 0.1—4.5 d/10a, and the variation was significant at Chifeng, 
Changzhi and Tacheng stations (P<0.05; Table 6). This result suggests that the observed changes 
in maize phenology, to some extent, were affected by SD adjustment. 

In M, experiment (Table 6), climate warming caused AD of maize to advance by 0.2-5.5 d/10a. 
Similarly, MD advanced by 0.6—11.1 d/10a at 5 stations. Due to climate-driven changes in maize 
phenology, the durations of VGP, RGP and WGP shortened by 0.2-5.5, 0.4-5.6 and 0.6-11.1 
d/10a, respectively (Table 6). 

Moreover, variation in SD had certain impact on maize phenology. The variation obtained from 
M>—M, reflected the impact of SD adjustment on maize phenology (Table 6). As shown in Table 6, 
the advancement of SD led to early AD and MD. Due to SD adjustment, AD occurred early by 
0.5-2.6 d/10a, and the variation was significant at Changzhi, Wuwei and Tacheng stations. 
Similarly, MD significantly advanced by 0.1-3.4 d/10a. SD adjustment caused VGP of maize to 
prolong at 4 of 5 stations (except for Fengning station), and duration of RGP slightly shortened by 
0.1-1.3 d/10a (Table 6). As a result, with advancement of SD of maize, duration of WGP 
lengthened at Chifeng, Changzhi and Tacheng stations, but it slightly shortened at Fengning and 
Wuwei stations (Table 6). 

From the comparisons of Mop—M2, cultivar shift during 1981-2010 significantly delayed 
anthesis of maize by 4.1, 3.4 and 3.8 d/10a at Fengning, Changzhi and Wuwei stations, 
respectively (Table 6). Due to cultivar shift, MD significantly delayed by 4.9-12.2 d/10a at 5 
stations (P<0.05). Furthermore, cultivar shift mainly prolonged duration of three growth periods 
(i.e., VGP, RGP and WGP) of maize at 5 stations. As a result, the durations of VGP, RGP and 
WGP lengthened by 0.24.1, 1.6-8.4 and 4.3—11.8 d/10a, respectively (Table 6). 

At 5 stations, climate warming in the past 30 years advanced AD and MD of maize by 2.4 and 
4.7 d/10a, respectively. Due to the impacts of climate warming on maize phenology, the durations 
of three growth periods (VGP, RGP and WGP) shortened by 2.4, 2.3 and 4.7 d/10a, respectively 
(Fig. 5). However, the negative impact of climate warming can be mitigated by the changes in 
crop management, such as SD adjustment and cultivar shift. In the past 30 years, SD of maize 
advanced in the study area. At 5 stations, the advance of SD caused AD and MD of maize to be 
advanced by 1.5 and 1.9 d/10a, respectively. However, SD adjustment caused VGP and WGP of 
maize to be prolonged by 1.2 and 0.6 d/10a, respectively (Fig. 5). Moreover, the negative impacts 
of climate warming on maize phenology were compensated as a result of introducing new maize 
cultivars with higher AT (thermal time) requirements. As shown in Figure 5, cultivar shift caused 
the AD and MD of maize to delay by 2.4 and 8.0 d/10a, respectively. Therefore, as a result of the 
impact of cultivar shift on maize phenology at 5 stations, the durations of VGP, RGP and WGP 
prolonged by 2.4, 5.8 and 7.1 d/10a, respectively (Fig. 5). 

Table 6 Observed and simulated variations in sowing date (SD), anthesis date (AD), maturity date (MD) and 


durations of three growth periods (VGP, RGP and WGP) of maize at 5 stations in northern China during 


1981-2010 
Phenological period Variable (d/10a) Chifeng Fengning Changzhi Wuwei Tacheng 
SD Mean value (DOY) 119.0 113.0 118.0 101.0 122.0 
Variation in Mo» -—3.0" —0.1 4.0" -1.7 —4.5* 
AD Variation in Mo -0.3 25 -0.4 -3.2™ -5.7" 
Variation in M; -0.2 -1.1 -1.8" -5.5™* -3.2™" 
Variation in M2 -1.1 -1.6 3.8" -7.0" 5.8" 
Variation in M)-M, -0.9 -0.5 -2.0° -1.5° -2.6" 
Variation in Mj,-M> 0.8 4.1” 3.4” 3.8™ 0.1 
MD Variation in Mop 4.4" 5.0" -0.5 -1.1 -0.7 
Variation in M, -0.6 -3.7"* -2.5* -11.1" -5.5™* 
Variation in M2 -1.7 3.8" 5.4" -13.3™ 8.9" 
= Variation in M)-M, -1.1 —0.1 -2.9™ —2.2" 3.4" 
Variation in Mj,-M> 6.1” 8.8” 4.9™ 12.2" 8.2" 
N VGP Variation in Mo 2.7 2.6" 3.5” -1.5 -1.1 
y Variation in Mı -0.2 z1 -1.8° 5.5" 39" 
Variation in M2 1.9" -1.5 0.1 4,9" -1.3 
Variation in Mo-M, 2.1" -0.4 1.9" 0.6 1.9" 
Variation in Mj,-M> 0.8 4.1 3.4” 3.4 0.2 
RGP Variation in Mo» 4.8” 2.5* 0.0 2.1 5.0” 
e? Variation in My -0.4 g” -0.7 -5.6" re 
pam Variation in M3 -0.6 27" -1.6 6.3" 3.3" 
NI Variation in M>-M, -0.2 -0.1 -0.9 -1.3 -1.0 
: Variation in M.,-M> 5.4" 5.2" 1.6 8.4" 8.3 
a WGP Variation in Mop 74° 5.17 35 0.6 3.9" 
4 Variation in Mı -0.6 aa as a" -5.5" 
> Variation in M2 1.3 4.5” -1.4 —11.2" 4.5" 
= Variation in MyM, 1.9 -0.8 1.1 -0.1 1.0 
Variation in Ma~ M2 6.1” 4.3™ 4.9* 11.8%” 8.4" 


Note: * and ** denote significance at P<0.05 and P<0.01 levels, respectively. VGP, vegetative growth period; RGP, reproductive growth 
period; WGP, whole growth period. DOY, day of year. 


MB Climate warming 
E] Sowing date adjustment 
L] Cultivar shift 


Variation of maize phenology (d) 
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Fig. 5 Effects of climate warming, sowing date adjustment and cultivar shift on maize phenology under 
different phenological periods at 5 stations in northern China during 1981-2010. Bar indicates standard error. 
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4 Discussion 


Under combined influence of climate warming, cultivar shift and SD adjustment, maize 
phenology experienced significant variation in NC during the past several decades. According to 
the results of this study, we found that variations of AD and MD of maize were opposite. When 
AD of maize slightly advanced, MD slightly delayed. The advancement or delay of maize 
phenology led to corresponding changes in the durations of different growth periods. Overall, 
three growth periods (i.e., VGP, RGP and WGP) prolonged in NC during 1981—2010. Generally, 
temperature is the most important climatic factor that influences crop phenology, because studies 
had shown that temperature impacts the growth of crop meristem that is the most important 
structure of crop growth (Hou et al., 2014). To some extent, the growth and developmental 
process of crop will be accelerated under climate warming (Qin et al., 2018). Therefore, 
increasing temperature often leads to earlier phenological period and a shorter growth period 
(Wang et al., 2015; Xiao et al., 2016a, b). In this study, climate warming during maize growth 
periods showed a significant increasing trend, which caused maize phenology (i.e., AD and MD) 
to advance, and eventually shorten the length of growth periods (including VGP, RGP and WGP). 

Generally, phenological periods of crop are jointly affected by changes in climate (mainly 
temperature) and improvement of agronomic management practice (including SD adjustment and 
cultivar shift) (Xiao et al., 2016b). Thus, it is unreasonable to simply attribute changes in 
phenology to climate change since agronomic management practices have also been changed (Liu 
et al., 2013). In this study, by using APSIM model, we evaluated the impacts of cultivar shift and 
SD adjustment on maize phenology in NC. During 1981—2010, cultivars were frequently shifted 
and new cultivars had a higher AT requirement. Overall, due to contribution of cultivar shift on 
maize phenology, durations of three growth periods (i.e., VGP, RGP and WGP) of maize 
prolonged. Although climate warming accelerates maize growth and shortens the growing periods 
of maize, cultivar shift significantly prolonged the growth period of maize in NC. As maize yield 
was determined by the length of growth period (Tao et al., 2014b), the cultivar shift with a high 
AT requirement, to some extent, mitigated the negative impact of climate warming. Therefore, 
introducing new cultivars with a longer growing period has been recommended as an adaptive 
measure to climate warming in the past several decades (Liu et al., 2013; Zhao et al., 2015; Xiao 
et al., 2016a, b; Wang et al., 2016; Liu et al., 2018). 

Optimizing SD is another approach for mitigating the adverse effects of climate warming. Thus, 
it is important to make appropriate SD adjustment for coping with climate warming. Generally, 
SD of crop mainly depended on climate conditions and soil moisture (Xiao et al., 2016b). Climate 
warming during past decades, to some extent, allowed for an advancement of SD of maize in NC. 
With the SD advancement, AD and MD of maize occurred early in NC. Studied had showed that 
SD adjustment exerted important effects on the length of growth period of maize (Lv et al., 2019). 
SD adjustment caused VGP of maize to prolong, while RGP was slightly shortened in this study. 
The above results are similar to the findings by Lv et al. (2019). They found that postponing SD 
of maize would extend the grain-filling stage (i.e., RGP). 


5 Conclusions 


This study analyzed the impact of climate warming, SD adjustment and cultivar shift on maize 
phenology in NC. The results showed that growth period of maize reduced due to the 
advancement of phenological periods, confirming that the dates of maize phenology shifted 
markedly in response to ongoing climate warming. However, this negative impact on maize 
phenology can be reduced by adopting new cultivars with higher AT requirements. Therefore, 
new maize cultivars with higher AT requirements and higher temperature tolerant should be used 
in response to climate warming. Moreover, it is needed to evaluate the relationship between crop 
phenology and yield, and further to adopt strategies to ensure food production under the 
background of climate warming. 
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